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Abstract 

In the collinear factorization of the form factor for the transition 7*7r° — > 7 the hard part 
contains double log terms as In 2 x with x as the momentum fraction of partons from to 1. A 
simple exponentiation for resummation leads to divergent results. We study the resummation 
of these In 2 x terms. We show that the In 2 x terms come partly from the light-cone wave 
function(LCWF) and partly from the form factor. We introduce a jet factor to factorize the 
In 2 x term in the form factor. To handel the In 2 x terms from the LCWF we introduce a 
nonstandard light-cone wave function(NLCWF) with the gauge links off the light-cone direction. 
An interesting relation between two wave function is found. With the introduced NLCWF and 
the jet factor we can re-factorize the form factor and obtain a new hard part which does not 
contain terms with In 2 x. Beside the renormalization scale /i the introduce NLCWF and jet 
factor have extra scales to characterize their x-behaviors. Using the evolutions of the extra 
scales and the relation we can do the resummation perturbatively in sense that the LCWF is 
the only nonpertubative object in the resumed formula. Our results with some models of LCWF 
show that there is a significant difference between numerical predictions with the resummation 
and that without the resummation, and the resummed predictions can describe the experimental 
data. 
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1. Introduction 



Predictions with perturbative QCD for an exclusive process can be made if it contains short- 
distance effects beside long-distance effects. In order to use the perturbative theory of QCD one 
needs to separate or factorize long-distance- and short-distance effects. Only the latter, which are 
characterized by a large energy scale denoted generically as Q, can be studied with perturbative 
QCD. It has been proposed long time ago that such a process can be studied by an expansion of 
the amplitude in 1/Q, corresponding to an expansion of QCD operators in twist[TJ[2]. The leading 
term can be factorized as a convolution of a hard part and light-cone wave functions of hadrons. 
The light-cone wave functions are defined with QCD operators, and the hard part describes hard 
scattering of partons at short distances. The hard part can be safely calculated with perturbative 
QCD in the sense that it does not contain any infrared(I.R.)- and collinear divergences. This is the 
so-called collinear factorization. In this factorization the transverse momenta of partons in parent 
hadrons are also expanded in the hard scattering part and they are neglected at leading twist. 

Although the hard part does not contain any of I.R- and collinear singularities, but as an 
perturbative expansion it contains large logarithms at higher orders of a s . These large logarithms 
are dangerous and can spoil the expansion in the sense that the expansion does not convergent. 
A resummation of large logarithms is often needed to have a reliable prediction. In this paper we 
study the resummation in the process 7*71"° — > 7. Theoretically, the process has been studied with 
QCD factorization extensively [3j IH HJ EJ El EJ [9]. Experimentally it has been studied too[lO]. With 
the collinear factorization the form factor F(Q) characterizing the process can be written as 

F(Q)~<£®i?{l + ^j. (1) 

In the above 4> is the light-cone wave function(LCWF) of 7r°, H is the hard part, and Q 2 is the 
virtuality of the virtual photon. The correction to the factorized form is power-suppressed and is 
proportional to A 2 /Q 2 . A is any possible soft scale characterizing nonperturbative scale, like Aqcd, 
the mass of ir etc.. This scale is about several hundreds MeV. In [8J the hard part in the collinear 
factorization has been calculated at one-loop, its two-loop result can be found in [9]. At one-loop 
level one finds that H contains a double log a s In 2 x, where x is the momentum fraction carried by 
a parton with < x < 1. It is expected that at n-loop level, it will contain a n In 2 ™ x. Those double 
log terms will become divergent when x is approaching to zero and can spoil the perturbative 
expansion of H. Although the double log terms are integrable with x in the convolution and give 
finite contributions, but they are significant corrections. The purpose of our study is to resum 
those double log terms. 

In order to resum those double log terms we need to understand their origin. Given the factorized 
form in the above, H will receive contributions from the form factor and the LCWF. When we use 
a finite quark mass to regularize the collinear singularities, we can show that the log terms come 
from the form factor and the LCWF. These terms can be re-factorized by introducing a jet factor 
and a nonstandard light-cone wave function(NLCWF), where the jet factor contains the double log 
from the form factor and the NLCWF doe not have double log. In introducing these two objects 
in the re-factorization, two extra scales, which will be explained later, are introduced. One of them 
is related to the In 2 x term from LCWF, while another is related to In 2 x in the form factor. The 
new hard part H from the re- factorization will not contain In 2 x. With evolution equations of these 
scales we are able to resum the In 2 x terms. Our approach is similar to the threshold resummation 
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in inclusive processes studied in [TT]. There exists an interesting relation between the two wave 
functions. In our approach, after the resummation of In 2 x, only the LCWF appears in the form 
factor as a nonperturbative object, other quantities can be calculated with perturbative QCD. With 
the knowledge of the LCWF we are able to give numerical results and to make a comparison with 
experiment. 

It should be noted that it is possible to resum the large log terms in exclusive processes by 
taking transverse momenta of partons into account [121 113j . By introducing transverse- momentum- 
dependent(TMD) light-cone wave-functions one can make TMD of kt factorization in terms of the 
TMD light-cone wave functions. However, the consistence of the factorization needs to be carefully 
checked beyond tree-level. Recently, this has been studied in a series of papers [TU [15], where a 
consistent definition of TMD light-cone wave functions is proposed and the TMD factorization is 
examined beyond tree-level for some simple cases. However, unlike LCWF's there is little knowledge 
about TMD light-cone wave functions. Therefore, it may be difficult to give detailed predictions 
and to compare with experiment. 

Our paper is organized as the following: In Sec. 2 we introduce our notations and explain the 
origin of In 2 x. In Sec. 3 we introduce NLCWF and present a one-loop study of the NLCWF. We 
also show that there is a perturbative relation between LCWF and NLCWF. In Sec. 4 we introduce 
our jet factor and derive the factorization formula with the jet factor and NLCWF. In Sec. 5 we 
give our resummation formula. We present our numerical results in Sec. 6, where a comparison with 
experiment is also given. Sec. 7 is our conclusion. 

2. Notations and the Origin of In 2 x 

We consider the process: 

tt° + 7* 7 (2) 

where 7r° carries the momentum P and the real photon momentum p. We use the light-cone 
coordinate system, in which a vector a M is expressed as a M = (a + ,a~,a±) = ((a + a 3 )/\/2, (a — 
a 3 )/\/2, or, a 2 ) and a\ = (a 1 ) 2 + (a 2 ) 2 . We take a frame in which the momentum P and p are: 

P^ = (P+,p- 0,0), p" = (0,p- 0,0). (3) 

We will consider the case that the virtual photon has the large negative virtuality q 2 = (P — p) 2 = 
—2P + p~ = — Q 2 . The process can be described by matrix element {jip, e*)| J£ m |tt (P)), which is 
parameterized with the form factor F(Q 2 ): 

<7(p,O|J£ m > (P)) = ie 2 e^e:P pPa F(Q 2 ). (4) 

e is the charge of proton, i.e., a = e 2 /(4-7r). 

In the collinear factorization the form factor can be factorized as 



1 + 



.A 2 , 
Q 2 



(5) 



where (j) is the LCWF of 7r°, H is a perturbative function or a hard part. Q u and Qd are the electric 
charge fraction of u and d quark in unit of e, respectively. <p is defined with QCD operators: 

= / ^e ifc+z ~(0|g(0)4(°°> 0h + l5 L n (oc,z-n)q(z-n)\7T°(P)), (6) 
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where the gauge link is defined along the light-cone direction = (0, 1,0,0) as: 



L n (oo, z) = P exp ^— ig s J dXn ■ G(Xn + z)j . (7) 
To find the hard part H, we replace the hadronic state with the partonic state: 



\n°(P))^\q(k q ),q(kg)), 



ki = ki 



m 



k2 



(k+, k-,o,o), 

(l-xo)P- 



kg 



kg 



= (k+, fc 3 -,o,o) 

x P+, 



(8) 



where we use a small but finite quark mass m to regularize collinear singularities. The form 
factor calculated with the partonic state will in general contain collinear singularities. The LCWF 
calculated with the partonic state will also have collinear singularities. If the collinear factorization 
holds, the singularities of the form factor and the LCWF will be the same so that the hard part H 
will not contain any collinear- and I.R. singularities. Here we examine this explicitly and show the 
origin of In 2 x. 



Figure 1: Feynman diagrams of tree- level contributions to the partonic scattering. The black dot 
denotes the insertion of the electric current operator corresponding to the virtual photon. 



With the partonic state, the LCWF at tree level reads: 

</> (0) (x,/x) = 5(x - x )(/) , 4> = v(k q )~f + ^u{k q )/P + . 



(9) 



At the leading order, the form factor receives contributions from diagrams in Fig.l. It is straight- 
forward to obtain the tree- level result from Fig.l: 



F(Q 2 



1 



la 



x 



1 



life 



3^2Q 2 (1 



3y/2Q- 



x 



(10) 



We will always use the notation u = 1 — u. Combining it with the tree-level result of the LCWF, 
we can obtain a factorized form for the form factor at tree-level: 



F(Q 2 



1 



3V2Q 2 



dxcj)(x, fj,) 



1 1 
- + - 

X X 



H^\x,Q^) = - + \. 
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At one-loop level, there are 12 Feynman diagrams, 6 of them are given in Fig. 2. The other 
6 diagrams are obtained from those in Fig. 2 by reversing the direction of the quark line, i.e., 
through charge conjugation. The diagrams in Fig. 2 represent the correction to Fig. la, and the 
other 6 diagrams for the correction to Fig. lb. Two corrections are related each other by charge 
conjugation. Hence we will need to study how the contributions from Fig. 2 can be factorized. 
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Figure 2: Feynman diagrams of the one-loop corrections to Fig. la. 



The contributions except that of Fig. 2a can be calculated in a straightforward way. We give 
gluons a small mass A to regularize I.R. singularities. The results are: 



F(Q 2 )W 



F(Q 2 )\2d 



F{Q 2 )U 



F{Q 2 )\2c 



F(Q 2 )\ 2f = F(Q 2 



la 



6tt 



In 



m 2 



2 In 



A 2 



F{Q 2 )\la- 



F(Q 2 )\la 



-a s 

3vr 
2a fi 



In ^— - In x + 1 

Q 2 



2 ln ^ +ln ^ + 2 lnX0 



1 



F(Q 2 )\la--^<- 

1 /i 2 m 2 

+ 2 W 



1 



Q 2 



7T 



-In xq + 2 lnxo lnxo + In — ^ hixo h 2Li2(xo) 



2 + xq 
2x 



In xq 



(12) 



These results can also be found in p3]. From the above, the In 2 x term comes only from Fig. 2c. 
The origin of this In 2 x is the following: The quark propagator connecting the vertex which emits 
the real photon carries the momentum (xP + , —p~, 0, 0) + 0(m 2 ). If x becomes small and goes to 
zero, the momentum becomes light-cone-like. If the momentum of the exchange gluon is in the 
region collinear to the — direction, after the loop integration a collinear singularity regularized 
by the small x appears. However, the region also overlaps with the infrared region where all 
components of the momentum are at order of xQ. This region generates an I.R. singularity which 
is also regularized by the small x. Therefore, the contribution from Fig. 2c contains In 2 x, one 
comes from the collinear singularity and another from the infrared singularity. We will show later 
that the dominant contribution containing In 2 x in these regions can be obtained by the eikonal 
approximation and can be factorized by using the method suggested in [16] . It has also been 
suggested by using a jet factor to absorb the In 2 x term [IT]. With a similar analysis one can show 
that the contribution from Fig. 2b does not contain the collinear singularity related to the — 
direction when x becomes small. It contains only the infrared singularity regularized by x. Hence 
it does not contain In 2 x as shown explicitly in Eq.(12). 

The one-loop correction of the LCWF is given by some of diagrams in Fig. 3 and Fig. 4. The 
diagrams with the gluon-exchange between gauge links give no contribution here because n 2 = 0. 
The one-loop part of LCWF is the sum: 



<jp-\x,n) = <j){x, /i)| 3o + 0(x, H)\zb + 4>{x, /^)|3c + 0(^, M)|4a + 0(^,M)|4c + 0(^,M)U(i + ^,M)|4/ ) (13) 
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Figure 3: Feynman diagrams of the one-loop corrections to LCWF. The double line stands for the 
gauge link 



The one-loop results can be found in p3]. They are: 

O'.s 



00, m) I 



la 



4>{x,n)\ 



2a, 



id 



6ir 



In At +21n 



m„ 



777 



{x^)\zc + (t>{x,n)\Af = — —(f) Q 9(x-xo) 

07T 



4>{x^)\?,b + 4>(x^)\ 



Ac 



A 2 



2as 
3vr 



<p 9(x - x) 



x {x - x 
x 

Xq(x - X ) 



m 2 (x - x ) 2 

m 2 (x-x f 
m 2~~2 



+ 



(14) 



Figure 4: Feynman diagrams of the one-loop corrections to LCWF. The double line stands for the 
gauge link 



With the above results the one-loop contribution of H can be determined as: 

HW(x ,Q,v) = FW(Q)- f 1 dx^ 1 \x,fi)H^(x,Q,fi). 

Jo 



(15) 



It is obvious that the contributions from Fig.2e and Fig.2f are already contained in the contribution 
of Fig. 4a and Fig. 4d of the LCWF, respectively. The contribution from Fig. 3a has a complicated 
expression. However, for the purpose of the factorization, we only need the contribution at the 
leading power of Q 2 . We have the result for the combination: 



Q 2 F(Q 2 )\ 2a 



2a s 1 
Sir x 



Inx 



Q 2 



l 



In — ^ — 1 H — In x 

/j z 2 



The convolution of other one-loop parts of LCWF reads: 

2a s 1 



Jo X 



+ 0(Q- 



i 1 



(16) 



3n XqXq 



In xq + 2 In xq In x + In xq In — ^ + 2Li 2 (xo) 



6 



7T 2 , II 2 

— + x hL—= + 2x 
3 m z 



dx 
o x 



(f>(x,fj,)\ 



3b+4c 



2a s 
3-kxq 



In 



(17) 



We note that the result in the first line can be used to subtract collinear singularities in Fig. 2c and 
that in the second line can be used for Fig. 2b. With these results one can extract the contributions 
of H from Fig. 2c and Fig. 2b: 



H^(x,Q,fi)\2 C = 
H^(x,Q,n)\ 2b = 




In 2 x — 2 In x In — (2 + x) In x 



Q 2 




ln- + lnx-4 . 



(18) 



Finally the one-loop part of H can be given: 



1 n 

H(x,Q, f i) = - + ^L- 

X 6TTX 



X O 2 
In 2 x - - lnx - 9 + In (3 + 2 lnx) 
x fl Z 



+ (x^x)+0(a 2 ) (19) 



From Eq.(17) we can see that the LCWF gives also a contribution with In 2 x to H. The origin 
of this double log is that we use the light-cone gauge link. With the light-cone gauge link the 
contribution from Fig. 3c and Fig.4f has a light-cone singularity beside a collinear singularity. The 
light-cone singularity is canceled in the sum. If we use gauge links with non light-cone vectors, the 
light-cone singularity will be regularized by the deviation of the vectors from the light-cone vector 
n. 

The obtained H behaves like xH ~ 1 + a s In 2 x/(37r) when x goes to zero. A resummation 
with a simple exponentiation does not work because of the +- sign in the front of the In 2 x term. 
Inspecting the one-loop part H one may chose /ias/i 2 = \fxQ 2 to kill the In 2 x term. However, for 
small enough x the scale becomes so small that perturbative QCD can not be used. It seems that 
one needs extra nonperturbative objects beside the LCWF to complete the resummation. We will 
show in our work that the resummation can be done without those extra nonperturbative objects. 

Before ending the section we would like to discuss the case if the dimensional regularization is 
used to regularize collinear singularities. In this case, the origin of In 2 x is different than that with 
a quark mass. However the hard part is the same and it is expected that quantities which are free 
from collinear singularities will not depend how the collinear singularities are regularized. 



3. Nonstandard Light Cone Wave Function 

As discussed before, one can use non-light cone gauge links to define nonstandard light cone 
wave functions. A possible definition is a straightforward generalization of Eq.(6): 

<Mx,C,/x)~ I ^ re lfc+ ^(0|g(0)Lt(oo,0) 7 + 75^(oo,z-n)g(z-n)|7r°(P)), (20) 
where the gauge link is 

L u (oo,z) = Pexp(-ig s J dXu ■ G(Xu + z)j , = (u + , u~ , 0, 0). (21) 
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This definition is gauge invariant. The defined NLCWF depends on an extra parameter 



u + u z 

We will take the limit u~ >> u + or Q — > oo. The limit C, — > oo should be understood as that we 
do not take the contributions proportional to any positive power of u + ju~ into account. It has no 
light-cone singularities as we will show through our one- loop result. 

At tree-level the NLCWF is the same as the LCWF. At one-loop level, there are contributions 
given by all diagrams given in Fig. 3 and Fig. 4. However, the contributions from interactions between 
gauge links will cause some problems, especially the contribution from Fig. 3d. It should be noted 
that the contributions from interactions between gauge links, i.e., those from Fig. 3d, Fig. 4b and 
Fig.3e, have no corresponding contributions in the form factor. A direct calculation shows that the 
contribution from Fig. 3d is not zero when x = or x = 1. When this contribution convoluted with 
the tree-level hard part it will lead to divergences. Therefore these contributions need to be 
subtracted and a modification of the above definition is needed. Since only interactions between 
gauge links are involved in these contributions, one can consider to use products of gauge links to 
subtract them. 




Figure 5: The one-loop contribution to products S of gauge links. 




Figure 6: The one-loop contribution to products S of gauge links. 
We consider the following products of gauge links: 

S(z~,u,v) = -^-Tr(0|Lj(0,-oo)L^(oo,0)L u (oo,2;~n)L t ,(z _ n,-oo)|0), 

* 'c 

S(z~,u,n) = j^Tr(0\Ll(0,-oo)Ll(oo,0)L u (oo,z~n)L n (z~n,-oo)\0), 

S{z~,n,v) = — Tr(0|4(0,-oo)Lt(oo,0)L n (oo,2-n)L„(^-n,-oo)|0). (23) 
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The vector v is taken as = (v + , v , 0, 0) with v + » v . The fourier transformed S is: 

f dz~ 



S(q + ,u,v)= J -±- e ^ +z -S(z-,u,v). 



(24) 



At tree- level all S"s are 1 in the z~-space or 5(q + ) in the g + -space. At one- loop level, they receive 
corrections from Fig.5 and Fig.6. It is interesting to note that there are certain relations between 
contributions of the three gauge link products. E.g., the contribution from Fig.5b to S(q + ,u, v) is: 



„, , .4 , f d 4 k , , u-v 

S(q ,u,v)\ 5b = i - 9s j —Jik +q )- 1 - 



iO u-k-iQ k 2 -X 2 +iO' 
It is interesting to note that under the limit v + » v~ and u~ >> u + : 



u • v 



1 



(v ■ k — i0)(u ■ k — iO) n ■ k 



n • v 



n • u 



v • k — iO u • k — iO 



(25) 



(26) 



where the first term corresponds to the contribution from Fig. 5b to S(q + ,n, v), and the second 
term corresponds to that to S(q + ,u,n), Hence we have: 



S(q + ,u,v)\ 5b = S(q + ,n,v)\ 5b - S(q + ,u,n)\ 5b . 



(27) 



The same result also holds for Fig.5c and Fig.6e and Fig.6f. With this observation we define the 
soft factor: 



S(z ,(u) 
S(q + ,(u) 



1 

2 

P 



[l + S(z ,u,v) — S(z ,n,v) + S(z , u, n)] , 
dz~ , 



I 



2vr 



(28) 



With above results we have up to one-loop level: 



2rv 

S((x-x )P + ,Cu) = 1 + — ^ 

07T 



1 



X - X J + 



[6(x - x) - 9{x - xq)] 



--S(x-x ) 



In 



- xo) 2 , ^ Clxl 



+ ln 



+ 0{a 2 s ). (29) 



The soft factor S only receives contributions from Fig.5d and from the self- interaction of a gauge 
link, i.e., from Fig. 6a, Fig. 6b, Fig.6d and Fig.6e. These contributions are the same of those from 
Fig.3d, Fig.4b and Fig.4c to the NLCWF, respectively. Therefore we can use this fact to subtract 
the contributions to 4> + from Fig.3d, Fig.4b and Fig.4e. It should be noted that the soft factor 
used here may be not unique. This non-uniqueness may be fixed through a study of higher orders 
and will not affect our one-loop results in this work. Also in the soft factor S one can take the 
limit v — > I so that S does not depend on v. 
We modify the definition of the NLCWF as: 



dz~ tk+z - <0|9(0)Lt(oo, 0) 7 +75^(oo, z~ n)q{z- n)\^ {P)) 



2vr 



(30) 



9 



Its individual one-loop contributions are 



2a, 



<f>+(x,0hb + 4>+(x,0W = "3^0 { 0(XO - X) ( - — 



<^ + (£,C)|3c + <Mz,C)|4/ = 



2a 



X 



Xq X — Xq 



In 



+(5(x - x ) 



1, H Z TT " 

2 ln «-T + 1 



^2 2 
> 00, 



—0o < 0(x - x ) ( — • 



X 



+(5(x - x ) 



Xo X — Xo 
2 



In 



c 2 4 



mi 



2 m C 2 x 2 



7T 

ir + 1 



(31) 



it should be noted that there is no term like (ln(x — xo)/(x — xo))+ in comparison with <fi + (x,n), 
hence it will not lead to any term with ln 2 x when convoluted with H^°\ The one loop result for 
d> reads: 



$\x,(,n) = =^{e(x -x)[-- 



2a 



3vr 



+S(x - Xq) 



X 



1 



In 



1 in " 



x x - x / + 

2 1 ,,2 



(■2 2 

^ + ^(x-xo)(^-- 



2 - C 2 x 2 



Mn ^ 



m 

7T 



1 



2 _ C 2 ^ 2 



+ 2 







Xo X — Xo 

■0 + Fig. la 



In 



mi 



+-^-S(x - x ) 

07T 



_ ln i^ +2 ln-|-4 



(32) 



the last line is from external- leg corrections. 

There is an interesting relation between LCWF and NLCWF. It reads: 



-(x, (,(*)=[ dyC(x,y,(,n)(f)(y,n), 
Jo 



(33) 



where the function C can be calculated with perturbative QCD and does not contain any soft 
divergence. From our results we have: 



C(x,y,C,b,fj.) 



S(x -y) + 
-6(y - x) 
+S(x - y) 



2a 



3vr [ 
1 



9{x - y) 



1 



' x ^ C 2 {x-y) r 



x-y \y 

2/rr. „.\2" 



I 1 



x ln C (x - y) 



x-y \y 



1, A* 2 1, A* 
— In n n H — In n „ 

2 CV 2 CV 



7T 



+ 0(a 2 ). 



With the function we define another function which will be useful later: 

C(x,(,fJ>) f 1 dy 



x 



f —C(y,x,£,n). 

Jo y 



(34) 



(35) 



which is just the convolution of C with the tree-level hard part from Fig. la. Again the function 
has a perturbative expansion: 



C(x,C,/x) 



1 - 



2a s (/x) 1 



3vr 



x 



-ln 2 x + lnxln ^2 + 2Li 2 (x) - 



7T 
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2 1 " C 2 .i- 2 



+ - In 

2 



(36) 



4. The Jet Factor and Re-Factorization 



After having studied the double log in LCWF, we need now to study how to factorize the double 
log In 2 xq from the form factor from Fig. 2c. As discussed before, the double log comes from the 
loop-momentum region where all components of the momentum carried by the gluon are at order 
of xqQ. One can use the eikonal approximation to expand the contribution in xq before the loop 
integration to obtain the dominant contribution. After some algebra we have: 



<7(P,Ol^mJ^,%>|2c 



4f7 2 f d 4 k 1 



/ 



P 



3 J (2vr) 4 fc 2 + i0 (k + k q - p) 2 - m 2 + iO -2k~ + iO 
(l(p,e*)\^. m .\k q ,k q )\ la + O(x° ), 



(37) 



where we omitted irrelevant factors. The eikonal propagator l/(— 2k~ + iO) comes from the quark 
propagator from the anti-quark after emitting the gluon. This suggests that we can replace the 
anti-quark line with a suitable gauge link. However if we take the gauge link along the light-cone 
direction, it will produce a light-cone singularity in the integration over k~ as indicated above. To 
avoid this we can take the gauge link along non light-cone direction. 

For our purpose we consider the following time-ordered product of gauge links with quark fields: 



S q (z) = l 1 ft{ 7 -<0|r[K t (0,-oo) g (0)g(z)^ i (z,-oo)]|0)} ) 
SM = J ■ 



(38) 



with = (xP + , —p~,0, 0). Without the gauge links, it is just a quark propagator. We first fix the 
vector v with v + » v~ as discussed above. The gauge link with u is needed to make S q (q) gauge 
invariant. The direction of u will be give later. At tree-level we have 



il (7 • q + rn) 



7 2 — m 2 



1 



(39) 



At one- loop there are corrections from diagrams given in Fig. 7. The dominant contribution 
from Fig.7c is proportional to the factor in [• • •] in Eq.(38) if the eikonal propagator l/(— IkT + iO) 
is replaced with l/(—2v ■ k + iO). Hence it will produce the same In 2 x as that in the form factor 
from Fig. 2c. However, there are contributions involved with interactions only between gauge links. 
They are those from Fig. 7d. Fig 7e and Fig.7f. If we take the direction u as = (0, 0, u 1 , u 2 ), the 
contribution from Fig.7d can be eliminated because of u ■ v = 0. The other two can be subtracted 
with expectation value of gauge links. We define the following jet factor as: 



q + f 

J(x,(-,,Q,n) = — J d 4 ze 



■iq-z 



Tr{ 7 -(0|T 


V,t(0,-oo)g(0) q{z)Vu(z 


-00) 


|0>} 


Tr(0|T 


Vj(0,-oo)Va(0,-oo)" 


|0> 



f = (xP + ,-p-,0,0), Q 2 = 2P + P ~, d = 



2v+(p- 



(40) 
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Figure 7: One-loop correction of S q . 



with the denominator the contribution from Fig.7e and Fig.7f are subtracted. The tree- level con- 
tribution to J is 1. The one- loop contributions are then from Fig. 7a, Fig. 7b and Fig. 7c. They 
are 



J(x,C 7 ,Q,/U)|7a 
J{x,(j,Q,V)\7b 

J(x,C^,Q,fi)\ 7c 



3^r 



In ^ 

hi 



+ 2 



3vr 



In 



c 2 

+ 2 + In > 7 



C 2 

In 2 ^ 7 



xQ 2 



vr 2 -4 



+ o(C 



7 ' 



(41) 



At the order we consider J does not depend on u 2 . We note that the In 2 x term from Fig. 7c is 
exactly that from Fig. 2c contributing to the form factor as expected. With the above result one 
can derive the following evolution equations which will be useful for our resummation: 







d In fi 2 

d 



J(x,C 7 ,<2,yu) 



3vr 



J(x,("f,Q,fJ>), 

1 -21n-L + l j J(x,C 7 ,Q,/x). 



With NLCWF and the jet factor we can write a factorized form for the form factor: 



HQ) 



dx 



^<p + (x,C,^)J( x X'y,Q,fJ')H(x,C,C'y,Q,^) +(x->l-x) 



(42) 



(43) 



and H (x) does not contain In 2 x explicitly. The leading term H is one in the above. With our 
NLCWF we have the convolutions corresponding to that with LCWF in Eq.(17), subtracted with 
the corresponding contributions from the form factor: 



Q 2 F(Q 2 )\ 2c 



Q 2 F(Q 2 )\ 2b 



— 4>(x,()\3c+4f 
X 



— Hx,C)\3b+4c 
X 



2 °s A 






onxo 


x V 


2 




iln 
2 


( 2 x 




Q 2 



Q 2 

— I --In 2 x + lnx ln— 2- + 2Li 2 (x ) 



C 2 



7T 



~2xV lnX0+ 6 



2a. 
37ra;o 



i, Q 2 

— in — = — 

2 C 2 io 



2^-1 

6 



(44) 
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it is clearly that all collinear singularities related to the quark mass are factorized into the NLCWF. 
With the jet factor we have: 



H(x,(,(j,Q,ij,) 



1 + 



37T I X 



In x In — r- + - In — 5- + In x In — r 

C 2 2 C- 



In + In H — In ^ 



c 2 



c 2 



+ 



1 



7 

X7T 2 \ 7T^ 

— + 2Li 2 (x)j -2+y-lnx 



3x , x < 
— In x In' 



2x 



2x 



xQ 2 

c 2 



Q 



1 



-lnx I In -Zj - 1 + -lnx I ^ + O(af), 



(45) 



as expected, for fixed (, C7 and \x there are no terms like ln 2 x. Also there are no terms like lnx 
without involving other log's. 



5. Resummation 



If we take ^ + as a nonperturbative object entirely, the resummation is really simple, in which 
we chose scales like /j,, ( and C7 so that there is no large log's. E.g., we can take those scales and 
obtain H which contains no large log's: 



f = Q2 = C 2 = C 2 



H(x,Q,Q,Q,Q) 



1 + 



X7T 



3vr 



+ 2Li 2 (x) 



7T 

2 + y-lnx 



3x x 

In x 

2x 2x 



ln 2 x- I lnx ^-1 + 2 ln ^)} + 



(46) 



and for the jet factor we use the evolution equation of £ 7 to express J at £7 = Q with that at 

C 7 = V%Q : 



J(x,Q,Q,Q) = exp 



a s (Q) 
3vr 



In 2 x + In x 



J(x, v^Q, Q, Q), 



J(x, v^Q, <2, Q) 



1 + 



3^ 



— In x — 1 — 7r z 



(47) 



it should be noted that J(x, y/xQ, Q, Q) has no In 2 x term. Only the single log lnx remains. Then 
for the form factor we have: 



F(Q) ~ / dx 
Jo 

+(x 1 - x)] . 



-<^>+(x, Q)J(x, \/xQ, Q, Q)H{x, Q, Q, Q, Q) exp J_ fHfl^) 
x I Sir 



In x + In x 



(48) 



Taking H J as a perturbative function, it does not contain In 2 x. The terms with In 2 x are re- 
summedin the exponential, but one needs the information of (j) to make predictions. 

With our results one can indeed resum In 2 x in the factorization formula with LCWF. One can 
use the relation between (j) and <fi to write another factorization formula for the form factor: 



F(Q) ~ /' 
Jo 



dx 



x 



4>(x, fi)C(x, (, n)J (x, C 7 , Q, (i)H(x, C, C 7 ; Q,f*) + (x ->l-x) 



(49) 
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and take C, J and H as perturbative functions. H is the same as H at one-loop level. If we expand 
the product CJH, we return to the standard collinear factorization discussed in Sec. 2. With the 
product its each part has a clear meaning. The ^-dependence in C will control the behavior of 
x — > in LCWF, while the Cy dependence in J controls that in the form factor. The evolution 
equations of C can be obtained from results in Sec. 3. They are: 



dln/i 2 



(50) 



For the resummation we first chose a scale [i\ in the factorization formula and use the /i-evolution 
to express C(x, C, fj,i) with C(x, (, fx): 



C(x,C,m) = exp<! -^|~hi ° s ^\ f 1 -^- + 1 H C(x,(,fi), 



3f3 a s (fii) V x 
where we have used the one-loop a s -running: 



(51) 



(52) 



We then use the (-evolution to express C(x, (, n) with C(x, Co, A*): 

C(x, C, /*) = exp | ^ + i) i n 1 1 c(x, Co, /*), 



(53) 



now we take Co = A* 2 / 21 , s0 that C(x, Co, A 4 ) does not have any log's: 



C(x,n/y/x,ll) = 1 



Q s (/i) 
3vr 



1 / vr 2 



1 vr 2 

- In x -\ h 2 

2 3 



+ 0(a 2 ). 



(54) 



With these steps, all log terms in C are resumed: 
C(x,(,(J-i) = exp ^ - 



8 q a (/i) _ 2a s (/i) xi^ 



3/3 a s (/ii; 



3vr 



A^ 



^ + 1) KV./'/v^-/'). (-^) 



Now we have the freedom to chose /x so that the exponential does not go to oo when x goes to 0. 
We can take \i fixed by: 



A* 



2 1 // 2 



a s (//) = xa s (//i), ln X2 = x ln A2' 



(56) 



and 



C(x,C,Mi) = exp j- — 
For a;-tOwe have now: 



lnx — ^-xa s (u,i) In x — x s + 1 
4vr a s (Q 



C(x,(,m) ~ exp 



3/3 



In x 



^. + l)\C(.r., l /^.f l ). (.77) 



(58) 
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and it goes to zero fast than any positive power of x. 
By taking the scales 

M? = Q 2 = C 2 -- 



(59) 



and using the (^-evolution to express J at £~ = Q 2 with J at ( 2 = xQ 2 , we obtain our resummed- 
form for the form factor: 



F{Q) 
S(x,Q) 



dx 



+(x 

8 
Wo 



-Mx, Q)J(x, \fxQ, Q, Q)H(x, Q, Q, Q, Q)C(x, [ij\fx, \x) exp {—Six, Q)} 
x 

1-x)], 



lnx — — xaJQ) lnx — x + 1 

4-7T 



+ 1 + 



.r 



« S (Q) r 



3?r 



In 2 x + In x 



8 | a fl (Q) 
3/3q 37r 



In x, for x — > 0, 



(60) 



in the above the product CHJ does not contain any log's, except J has a single log lnx. All other 
logs, like ln 2 x,etc, are resummedin S. Since we only used one-loop evolutions, for consistence we 
should neglect higher orders in a s in the product. Therefore, we have the one-loop resummedform 
factor: 



HQ) 



i 



3V2Q 2 
3Q 2 



dx 



-4>(x, Q) exp {— S(x, Q)} + (x — ► 1 — x) 



f 1 

/ dx— 4>(x, Q) exp {— S(x, Q)} . 
Jo x 



(61) 



This form can be used if one can get 4>(x, Q) easily at the large scale Q. If one only knows the wave 
function at lower scale, but does not want to solve the evolution equation for the wave function to 
get it at a higher scale, one can first evolute everything at a lower scale //o where the wave function 
is known or modeled, then to a higher scale [i\ and perform the resummation. We get in this case: 



\/2 f ^ 1 f 8 

F (Q) = ^TyI / dx-<t>(x,no)exp \ -S(x,Q) - — In 
3Q Z Jo x L 3po 



3Q 2 

where we have used: 



a s (Q) 

Q!s(Mo) 



lnx + 



djCHj) _ 2aM( lnx , 3 
d In /i 2 37r 



CHJ) . 



(62) 



(63) 



6. Numerical Results and Comparison with Experiment 

We will use our resummation formula in Eq.(62) and Eq.(63) to give our numerical results. In 
our formulas the nonperturbative input is the LCWF. The LCWF has the asymptotic form if fi 
goes to oo: 

<j>{x,n) = 6x(l-x)U + ---, (64) 

where • • • stand for terms which are zero in the limit \i — ► oo. The LCWF can be expanded with 
Gegenbauer polynomials [T]. A model for 4> has been proposed by truncating the expansion [18] : 

0(x,/z) = 6/^(1 - x) (l + </» 2 (/i)C 2 3/2 (2x - 1)) , (65) 
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where <fe(W)) is determined by QCD sum-rule method at /xq = lGeV[18j: 



02(^0 = iGeV) = 0.44. 



(66) 



We will use these two types of LCWF to give our numerical results. We will use Eq.(61) with 
the asymptotic form of 4> to make our numerical predictions. For LCWF given in Eq.(65) we use 
Eq.(62). We take the A-parameter as A = 237MeV. Our numerical results do not strongly depend 
on the value of A. There is only a little change if we change A from lOOMeV to 300MeV. 



0.35 




0.05 - 



Figure 8: Numerical results with experimental data. The curve A and curve As are obtained by 
using the asymptotic form without and with the resummation, respectively. The curve B and curve 
Bs are obtained by using the LCWF in Eq.(65) without and with the resummation, respectively. 
The experimental data are taken from the second reference in |10j . 

Our numerical results are given in Fig. 8, where experimental results from CLEO in [lOj are 
also given. From Fig. 8 we see that with the two types of LCWF the resummation has significant 
effects. The resummation can reduce the form factor predicted without the resummation at the 
level of 40% or more. Using the LCWF given in Eq.(65) with our resummation the experimental 
results can be well described for Q 2 > 3GeV 2 . 

7. Conclusion 

In the collinear factorization the form factor of the transition 7*71"° — > 7 can be written as a 
convolution of a hard part and LCWF. The hard part contains double log terms as In 2 x at one-loop 
level and is expected to have terms \v? n x at order of a". A resummation of these terms with a 
simple exponentiation can not be done because it results in divergent results. In this work we have 
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studied the resummation of these In 2 x terms. With a small but finite quark mass as the regulator 
of collinear singularities, we have found that the In 2 x terms come partly from the light-cone wave 
function and partly from the form factor, as discussed in Sec. 2. To handel these terms, we first 
introduce a nonstandard light-cone wave function with the gauge links off the light-cone direction. 
This introduces an extra scale in the NLCWF beside the renormalization scale /u. The deviation 
from the light-cone direction will regularize light-cone singularities in each contributions. This 
fact leads to that the NLCWF will not deliver any term with In 2 x to the hard part, if one uses 
the NLCWF to perform the factorization. As the next, we introduce a jet factor to factorize the 
In 2 x term in the form factor. The jet factor also contains an extra scale beside \i. This extra 
scale controls the x-behavior of the jet factor. Our re-factorized formula for the form factor is a 
convolution with the NLCWF, the jet factor and a hard part. The hard part does not contain 
terms with In 2 x. 

We have found that there is an interesting relation between the LCWF and the introduced 
NLCWF. The relation can be determined with perturbative QCD and is given at one-loop level 
in this work. With this relation we are able to show that the In 2 x can be resumed and the 
nonperturbative object in the resummed formula is only the LCWF. With the knowledge of LCWF's 
we are able to get numerical predictions. In performing the resummation of the double log we have 
used the concept of QCD factorization and worked out every quantity explicitly at one-loop level. 
It is possible to extend our work to the resummation of the remaining single log terms and beyond 
one-loop level. 

Our numerical results show that the effect of the resummation is significant. There is a difference 
at the level of 40% or more between the predicted form factors with and without the resummation. 
In comparison with experiment we find that the numerical predictions by using the LCWF in 
Eq.(65) with the resummation are consistent with the experimental data. 
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